Kondo effect in optics of quantum dots 
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Anderson impurity model for semiconductor quantum dot is extended to take into account both 
particle and hole branches of charge excitations. It is shown that in dots with even number of 
electrons where the Kondo effect is absent in the ground state, novel midgap excitonic states emerge 
in the energy spectrum due to Kondo-type shake-up processes. The relevance of the model to rare 
earth ions adsorbed on metallic surfaces is discussed. 



The Kondo effect was formulated initially for a spin 
1/2 impurity coupled to metallic environment. Later on, 
various realizations of this effect were studied in more 
complicated physical systems. Examples are the Kondo 
effect influenced by crystal field excitations, the multi- 
channel Kondo problem, electron scattering by two-level 
systems, etc [Q. Usually, the scattering center is repre- 
sented by its spin (pseudo-spin) degrees of freedom. In 
this work we consider yet another generalization of the 
Anderson impurity model, where the impurity possesses 
both electron and hole branches of single-particle excita- 
tions as well as two-particle (excitonic) states. Discrete 
excitonic states do not show up in a conventional sit- 
uation of an impurity immersed in a conduction band, 
even when such an impurity has internal charge degrees 
of freedom (e.g., rare-earth impurities with two unfilled 
electron shells Q). However, long-living excitons can ex- 
ist when the scattering center is spatially isolated from 
the conduction electrons, like in semiconductor quantum 
dots coupled by tunneling to metallic leads or in atoms 
adsorbed on metallic surfaces Q. It is well known ||] 
that the Kondo scattering could results in resonant tun- 
nel current. Recently the Kondo features were observed 
in the conductance of quantum dots || . We demonstrate 
in this work that Kondo-type shake-up process leads to 
an exciton formation in a quantum dot with even num- 
ber of electrons (where the Kondo effect is absent in the 
ground state). 

To reveal new qualitative features which appear in the 
Kondo problem due to excitons, we consider the following 
model of semiconductor quantum dot: 



H — Hd + Hi, + H t — Hq + H t . 
The first term, Hd, describes the isolated dot 



(1) 



H d = E c n c - E v n v + ^-n c (n c - 1) + -^-n v (n v - 1), 

(2) 

where a is the spin projection, n c = d\ a d ca is the 
occupation number for the electrons at an empty con- 
duction (c) level, n v = J2a™v<7 = dyadic is the occu- 
pation number for holes at a filled valence level E v < Ef- 



We define the gap as an excitation energy of the electron- 
hole pair, A = E c — E v . The band electrons in the leads 
with the Fermi energy ep and the tunneling between the 
leads and the dot with the tunneling amplitudes V^ c v ^ 
are described by the second and third terms in (0) re- 
spectively, 



ka j=c,v 



(3) 



No Kondo coupling is expected in the ground state 
and the strong Hubbard repulsion U CyV suppresses dou- 
ble occupation of electron or hole levels, so that only the 
states with singly charged dot 



Ifcc) = 



\kv) 
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4jv,\0), (4) 



are admixed by tunneling to the neutral state |0) (the 
filled Fermi sea and valence level E v ) of the isolated dot, 
and one can consider this admixture as a small first order 
perturbation (spin 1/2 is considered). There is also no 
Kondo effect in photoexcitation spectra 

D° + hv^D ± + e,h 

where D°,D ± are the ground state of the dot, and the 
charged states with a hole v and excess electron c, re- 
spectively; e(h) stands for the electron (hole) promoted 
to the leads due to absorption of the light quantum hv. 
We will show, however, that the Kondo- like processes de- 
velop in the spectrum of electron-hole excitations of the 
quantum dot given by the operator 



Mo) = s+ |o> . 



(5) 



These states can be excited by means of photon ab- 
sorption/emission processes which are described by the 
Hamiltonian H' = a Pij<d\ a dja exp(-iujt), where Py 
is the matrix element of the dipole operator P. 



1 



The optical line shape is given by the Kubo-Greenwood 
formula , 



Wihv) ~ lm-(i\PR(hu)P\i) 

7T 



(6) 



where R(z) = (z — H)^ 1 , and (i| . . . |i) means averag- 
ing of initial state over the equilibrium ensemble. Due to 
Hubbard exclusion mechanism only the states |0) survive 
in the brackets of response function (j^) , and one has to 
calculate the retarded excitonic Green function 



G? e (z) = -i dte izt 6(t)(\[B^t)B(0)]\) 



(7) 



Employing the operator identity R = R + R H t R for 
the resolvent R we construct the perturbation series in 
H t where the first term is R ee = (e|i?|e). Having in mind 
the fact that the strong coupling limit T -C Tjc hardly 
can be achieved (see below), we restrict ourselves by the 
weak coupling limit of T exceeding the Kondo tempera- 
ture Tk- 

In a simplest approximation one neglects the interme- 
diate states with multiple electron-hole pairs in the leads 
(see, e.g., 0). Then the closed system of equations for 
the matrix elements R~ e = (s\R~\e) where (s\ are the 
states admixed to the exciton by the tunneling interac- 
tion, is obtained. The structure of the system is illus- 
trated in fig.l, and its solution for the excitonic Green 
function has the form R~ e (e) = [z - E (e)][X>(e)] _1 . The 
poles of this function are determined by the equation 



V(z) = dct 
The self energies are given by 



z-H (z) -Eoe(z) 
-E 0e (z) z-A-£ e (z) 



= 0. 



(8) 



_ So = ^ 
^2£ 0e = E" 



1+ 



2E P = E, 



where 



E 



W*Wlf k 



E 



W*Wlf k 

z + E v - e k ' 



(9) 



(10) 



fk is the equilibrium distribution function for lead elec- 
trons, fk = 1 — fk- The tunneling matrix elements w are 
defined as (0\H t \kc) » V2(e\H t \kv) « w c , (0\H t \kv) w 
—\/2(e\Ht\kc) w*. Then integration gives 



RcE 



where A r 



E„ 



2tt n 



D 2 



(11) 



e F , A„ = e F - E v , Tji = npwjwi, p is 
the density of states at the Fermi surface, Ep is the Fermi 
energy and D is the conduction electron bandwidth. The 
signes ± correspond to c,v, respectively, z — e + is. 

As was mentioned above, the renormalization of the 
ground state at e ~ is a trivial second-order effect in 



the case of E c > ep, E v < ep (even number of electrons 
in the ground state of the dot). However, novel features 
appear in the excitation spectrum at e ~ A C ,A„. In 
a process of tunnel relaxation the electron and the hole 
with finite energy e induce Kondo-like peaks at the Fermi 
level. The physics of " Kondo-excitons" can be illustrated 
by considering two limiting cases, 
(i) Symmetric configuration w c = w v = w, A c = A„ 
A 2 In this case Eo e vanishes identically, and the secular 
equation becomes 



e- A-E e (e) = 



(12) 



It has a solution at e ~ A = A + E e (A), correspond- 
ing to the normal exciton. Besides, ReE e (e) diverges at 



T K = Dexp(-7rA/2F) 



(13) 



As a result a peak arises in ImG ee at e around A/2 which 
corresponds to maximum of light absorption/emission 
at this energies. When the particle-hole symmetry is 
slightly violated (S ^ 0), A c = A/2 — d, A v = A/2 + S, 
S <C A/2, this midgap peak disappears with encreasing S 
due to cancellation of singular terms with opposite sign 
in eq. (|[). We will see below that this state is also fragile 
against the lifetime effect. 

(ii) Strongly asymmetric configurations, A c -C A^, w c ^ 
w v . This case is closer to the situation in real systems 
where the electrons are usually less localized then the 
holes. In the extreme case the level E v is below the bot- 
tom of the band e k , (" Schrieffer- Wolff (SW) limit" @), 
whereas the electron is in resonance with the band elec- 
trons. The secular equation (0) may be rewritten as 



2(e-A) = E^ 



(e 



E-ut> ) 



(14) 



To find the resonance solution at e ~ A„, we neglect the 
smooth contributions E^-(e) in comparison with the sin- 
gular self energies E~ (e), and use approximate value of 
Y<~ V (A V ) « —r/Ac in the denominator of the ratio in r.h.s. 
(rj = (w v /w c ) 2 -C 1). Then we have 



2(e-A) wE" + V \J:-f(A v+V A c 



(15) 



Like eq. (|l^) , this equation has a Kondo-like pole at 
e = A V ,T = Tk, where Tk — D exp(— 2i:A c /T cc ) and 
T cc « r cc (l — ?]A C /A V ). Repeating the procedure for 
the resonance at e ~ A c , we leave in ( p^ ) only the 
terms E^ and find the midgap peak at e = A c , with 

T K = Dexp(-2nA v /T vv ), T vv w T vv nA c /A v . 

Thus, the Kondo- type processes manifest themselves 
as a shake-up effect with a shake-up energy of A VtC , 
i.e. as a final state interaction between the (e, h) pair 
in the dot and the Fermi continuum in the lead. The 
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T-dependent logarithmic singularity in excitonic self en- 
ergy is a precursor of " orthogonality catastrophe" in close 
analogy with the corresponding anomaly in a d-electron 
self energy in the conventional Anderson model 0. In 
the latter case the Kondo peak transforms to undamped 
Abrikosov-Suhl resonance in a ground state M. However, 
this is not the case for the Kondo exciton because of the 
finite lifetime 77 of the (e, h) pair. The most important 
contributions to t\ are given by the same tunneling pro- 
cesses which are responsible for the very existence of the 
midgap states. To take them into account one should 
include the states with (e, h) pairs in the leads in the 
Green function expansion. These states appear in 4th 
order of the perturbation theory. In a non-crossing ap- 
proximation (NCA) they result in renormahzation of the 
self energies f|fl0[ , 

" D f{e)de 



d e - E c 

D 



f(e)de 



d e 



(16) 



E v -e-BT x ( e ) 
where B^ are the integrals similar to (^) . In particular 



B. 



b: 



f(e)de' 



-D 
D 



2T r 



f(e)de' 



—D 



2r, 



■£' - e 



■£'-£■ 



A 



(17) 



In a symmetric case ( |l2| ) the self energy i? ± (e) has 
imaginary part w 2r at e ~ A/2. Thus, the Kondo pro- 
cesses initiated by one of the partners in the electron-hole 
pair are killed by the damping of its counterpart due to 
the same tunneling mechanism. More interesting is the 
asymmetric configuration described by eq. ([lj). In this 
case the singularity of Sj c in electron channel at e ~ A„ 
survives in a SW limit for a hole [E v below the bottom of 
conduction band, Im_B+(ei?) = 0], whereas the electron 
lifetime given by lmB~ v ~ r cc kills the hole peak like in 
the symmetric case. The electron midgap state survives 
also when E v has a finite width, provided < Tk, i.e. 

2^A c /r cc <i nJ D/r 

Thus, in the case (ii) the main peak of optical tran- 
sition at hv = A is accompanied by a satellite peak at 
hv A„. The form of this peak is determined by eq. 
(|), i.e., by 



—ImG ee (hv) = —Ira 

IT 7T 



hv — A — Yi ee (hv) + ih.T l 



(18) 

where ~E~ e is given by the r.h.s. of eq. (fl5|).The lineshape 
W{hv) strongly depends on T/f K and 77 [§0]. 

Two systems where this theory can be applied are sug- 
gested below. The first one is an ensemble of semicon- 
ductor quantum dots (e.g., the nanosize Si clusters em- 
bedded in amorphous SiC>2 matrix H). Luminescence of 



confined excitons was observed in these clusters || , and 
the tunneling current through this system exhibits the 
Coulomb blockade effect |^Jll|. In the Coulomb block- 
ade regime the energy U which enters Hd, is introduced 
as U = E^ n+1) -El n) +e 2 /C eff , where E^ is the energy 
level of the n-th quantum state of the electron (hole) in 
the empty well (n = 1 in our case), and C e // is an effec- 
tive capacitance of the barrier layers. Experimental |ll]] 
and theoretical jl2| estimates of C e / / for Si nanoclusters 
give U « O.It 0.3 eV, hence the general condition for 
realization of the Kondo effect, T -C \Ei\,Ef + U, is satis- 
fied in these systems. However, the experimental value of 
U is small compared with A ~ 1.2-5-1.3 eV pl|. Taking 
into account finite U means inclusion of doubly occupied 
states 1 2c), \2v) in the set (^), (||). In close analogy with 
the conventional Anderson model ||] one expects redistri- 
bution of the spectral weight of neutral states |0) and |e) 
in favor of the states |2c), \2v), and increase of Tk with 
decreasing U. However, the inequality Tk -C A ensures 
the existence of the midgap states. 

Another important condition for the appearance of 
midgap states is long enough excitonic lifetime r e ^> 
H/Tk- In addition to tunneling contribution T discussed 
above, one should also take into account the electron- 
hole recombination in the dot resulting in a width 7 of 
the excitonic level. This damping gives the contribution 
to imaginary part ~ r7/7rA <C T in eqs. (12) and ( |l4| ) for 
a Kondo pole. In any case, the experimentally estimated 
value is r e ~ 10~ 6 s for the singlet exciton |^|. Therefore 
the Kondo-type processes can survive in these systems if 
Tk ^> 10~ 9 eV, which is a realistic condition. We there- 
fore believe that our model qualitatively describes those 
properties of Si nanoclusters which play a crucial role in 
the formation of Kondo-type states in optical spectra. 

The second possible realization is mixed valent rare 
earth atoms adsorbed on a metallic surface. It is known 
that the Anderson model can be applied to adatoms with 
strongly interacting electrons (see [Q for a review). In 
this case H t (||) corresponds to covalent bonding be- 
tween an adsorbed atom and a substrate, V is the cor- 
responding hybridization integral between the electrons 
in adatom and those in the nearest sites of the metallic 
surface layer, U is the intra-atomic Coulomb repulsion 
which prevents charging of the adatom in the process of 
chemisorption. The model was originally proposed for 
hydrogen atoms adsorbed on surfaces of transition met- 
als |l^]. Later on, the possibility of Kondo-type spin 
polarization of substrate electrons around adatom spin 
in the case of U ^ V 2 / D was discussed M]. The most 
promising candidates from the point of view of excitonic 
effects are the adatoms with unstable valence, e.g., Sm, 
whose ground state electronic configuration is 4/ 6 6s 2 . 
The Sm atoms can be adsorbed on surfaces of transi- 
tion metals (Ni, Co, Cu, Mo). In the process of adsorp- 
tion, Sm loses its s-electrons and exists in two charged 
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states Sm 2+ and Sm 3+ depending on the concentration 
of Sm ions on the surface. In particular, the isolated ions 
Sm 2+ (4/ 6 ) are observed on Mo surface at low submono- 
layer coverage []l7j . The unfilled 4f shell forms a resonant 
f-state close to the Fermi level of the metal. The excited 
5d state forms another level above Ef- Thus we arrive 
at a two-level system described by our Hamiltonian Hd 
(||). Since the ground state term of the configuration 
4/ 6 is a singlet 7 Fq, one cannot expect the Kondo cou- 
pling for such adatom. However, in the course of virtual 
transitions between the adatom and the substrate the 
states \kv) and \kc) arise with excess electron above ef 
(configuration 4/ 5 efe) and a hole below Ep (configuration 
4/ 5 5d/ife). According to our calculations, one can excite 
not only the conventional atomic excited state with en- 
ergy A = E(4f 5 5d) - E(4f 6 ) but also the midgap states 
with energy close to A' = S(4/ 5 cf) — i?(4/ 6 ) where e_p 
stands for the electron on the Fermi level of the substrate. 

To summarize, this work suggests a generalization of 
the Anderson impurity model which takes into account 
excitonic degrees of freedom. The model exhibits Kondo 
effect in the excited state, despite its absence in the 
ground state. It is predicted that satellite excitonic peaks 
of a Kondo origin can be seen in the optical spectra 
of semiconductor quantum dots or rare earth atoms ad- 
sorbed on metal surfaces. 

We thank A. Polman and M. Brongersma for stimu- 
lating discussions of the optical properties of Si nano- 
clustcrs. This research is supported by Israeli Academy 
grants " Nonlinear Current Response of Multilevel Quan- 
tum Systems" and "Strongly Correlated Electron Sys- 
tems in Restricted Geometries" . 



[1] D.L. Cox and A. Zawadowski, Adv. Phys., 47, 599 
(1998). 

[2] K.A. Kikoin and S.B. Burkatovskii, Journ. Moscow 

Phys. Soc.3, 139 (1993). 
[3] Q. Ye et al, Phys. Rev. B44, 1806 (1991). 
[4] T.L. Einstein et al, in Theory of Chemisorption, edited 

by J.R. Smith (Springer, Berlin, 1980), p. 184. 
[5] L.I. Glazman and M.E. Raikh, JETP Lett. 47, 452, 

(1988); T.K. Ng and P.A. Lee, Phys. Rev. Lett., 61, 1768 

(1988). 

[6] D. Goldhaber-Gordon et al, Nature 391, 156 (1998); 

S.M. Cronenwett et al, Science, 281, 540 (1998). 
[7] C. Lacroix, Journ. Phys. Fll, 2389 (1981); N.S. 

Wingreen and Y. Meir, Phys. Rev. B49, 11040 (1994). 
[8] A.C. Hewson, The Kondo Problem to Heavy Fermions, 

Cambridge, University Press, 1993. 
[9] K.S. Min et al, Appl. Phys. Lett. 69, 2033 (1996). 
[10] P. Coleman, Phys. Rev. B29, 3035 (1984); N. Bickers, 

Rev. Mod. Phys. 59, 845 (1987). 
[11] E.N. Nicollian and R. Tsu, J. Appl. Phys. 74, 4020 

(1993). 

[12] D. Babic et al, Phys. Rev. B45, 14150 (1992). 
[13] D. Babic and R. Tsu, Superlattices and Microstruct. 22, 
581 (1997). 

[14] M. Brongersma, Ph. D. Thesis, Utrecht, 1998. 

[15] D.M. Newns, Phys. Rev. 178, 1123 (1969). 

[16] J.R. Schrieffer, in Dynamic Aspects of Surface Physics, 

Rend. Scuola Enrico Fermi, Varenna, v. 58 (1974). 
[17] A. Stenborg et al, Phys. Rev. B40, 5916 (1989). 



Figure Caption 

Fig. 1. Building blocks of the Green function expansion and the secular equation (||). The arrows indicate the 
tunneling processes which connect different states of this set. 
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